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Introduction 

One of the main algebraic structures related to 1 + 1-dimensional integrable PDE of the 
form 

Ut = F(u, u x ,...,u n ) , n>2, Ui = d l x u(x,t) (0.1) 
is an infinite hierarchy of commuting flows or, the same, (generalized) symmetries PP 

u u = Gi(u, ...,u mi ) • (0.2) 

We identify the symmetry ()0.2j) with its right hand side G{. The symmetry G satisfies the 
equation 

D t (Cf)-F m (G) = 0, (0.3) 
where D t stands for the derivation in virtue of (jO.lj) and F* denotes the Frechet derivative of 

The dual objects for symmetries are cosymmetries which satisfy the equation 

D t (g) + F t M = 0, 

where F* is the differential operator adjoint to F*. The product gG of any cosymmetry g and 
symmetry G is a total x- derivative. It is well known pQ that for any conserved density p the 
variational derivative is a cosymmetry. 

The simplest symmetry for any equation (jO.lj) is u x . The usual way to get other symmetries 
is to act to u x by a recursion operator 71. By definition, the recursion operator is a ratio of two 
differential operators that satisfies the identity 

[D t - F*, 11}= Tit- [F*, 71] = 0. (0.4) 

It follows from (jO.Hj) and (j0.4)l that for any symmetry G the expression 71(G) is a symmetry as 
well. 

Most of known recursion operators have the following special form 

k 

H = R + ^2G i D- 1 g i) (0.5) 

i=l 

where R is a differential operator, G^ and <7, are some fixed symmetries and cosymmetries 
common for all members of the hierarchy. For all known examples the cosymmetries <?, are 
variational derivatives of conserved densities. Applying such operator to any symmetry we get 
a local expression, (i.e. a function of finite number of variables u, u x , ...«,,...) since the product 
of any symmetry and cosymmetry belongs to Im D x . Moreover, a different choice of integration 
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constants gives rise to an additional linear combination of the symmetries G±, . . . ,Gk- Probably 
for the first time ansatz ()0.5|) was used in [2]. We call recursion operators (j(J.5|) quasilocal. 

Most of known integrable equations (jO.lJ) can be written in a Hamiltonian form 

-*$)■ 

where p is a conserved density and TC is a Hamiltonian operator. It is known that this operator 
satisfies the equation 

(D t -F*)H = H(D t + Fi), (0.6) 

which means that H, takes cosymmetries to symmetries. Besides (|0.6j) the Hamiltonian operator 
should satisfy certain identities (see for example, |T]) equivalent to the skew-symmetricity and 
the Jacobi identity for the corresponding Poisson bracket. It is easy to see that the ratio TL2H1 1 
of any two Hamiltonian operators is a recursion operator. 

As the rule, the Hamiltonian operators are local (i.e. differential) or quasilocal operators. 
The latter means that 

m 

H = H + Y,G i D- 1 G i , (0.7) 

i=l 

where H is a differential operator and Gi, Gi are fixed symmetries. It is clear that acting by 
the operator (j0.7)l on any cosymmetry, we get a local symmetry. 

For example, for the Korteweg-de Vries equation 

Ut = u xxx + 6 u u x 

the simplest recursion operator (see |3j) 

K = D 2 X + Au + 2u x D- x (0.8) 

is quasilocal with k — 1, G\ — 2u x , and g\ — 1 . This operator is the ratio of two local 
Hamiltonian operators 

Ux = D x , H 2 = Dl + AuD x + 2u x . 

For the systematic investigation of quasilocal (weakly nonlocal in terminology of jl]) operators 
related to the Korteweg-de Vries and the nonlinear Schrodinger equations see jU E] . 

It is possible to prove that for the Korteweg-de Vries equation the associative algebra A of 
all quasilocal recursion operators is generated by operator (|0.8|) . In other words, this algebra 
is isomorphic to the algebra of all polynomials in one variable. 

Our main observation is that it is not true for such integrable models as the Krichever- 
Novikov and the Landau-Lifshitz equations. These equations play a role of the universal models 
for the classes of KdV-type and NLS-type equations. It appears that other integrable equations 
from these classes are limiting cases or can be linked to these models by differential substitutions 

PEj. 



3 



It turns out that for these models, known to be elliptic, the algebra A is isomorphic to the 
coordinate ring of the elliptic curve. Since the algebra A is defined in terms of the equation 
(jO.lj) only, this is an invariant way to associate a proper algebraic curve with any integrable 
equation and, in particular, to give a rigorous intrinsic definition of elliptic models. 

This paper is organized as follows. In Section 1, we consider the Krichever-Novikov equation 
[Hj, which is the simplest known one-field elliptic model. We show that there exist two quasilocal 
recursion operators of orders 4 and 6 related by the equation of elliptic curve. These recur- 
sion operators are ratios of the corresponding quasilocal Hamiltonian operators. The simplest 
quasilocal Hamiltonian operator of order —1 

H = u x D~ x u x 

for the Krichever-Novikov equation was found in From the results of this paper it follows 
that this equation has also a Hamiltonian operator of order 3. In Section 1 we present one more 
quasilocal Hamiltonian operator of the fifth order for the hierarchy of the Krichever-Novikov 
equation. It seems to be interesting to investigate this multi-Hamiltonian structure for the 
Krichever-Novikov equation in frames of the bi-Hamiltonian approach [Oj. 

In Section 2 we obtain similar results for the Landau-Lifshitz equation. It turns out that 
for this model there exist quasilocal recursion operators of orders 2 and 3 related by an elliptic 
curve equation. The corresponding quasilocal Hamiltonian operators are also found. 

In Sections 1,2 we have used the direct way to construct recursion and Hamiltonian oper- 
ators based on the ansatzes (|().5|) and (|().7jl . On the other hand there are several schemes that 
use L — A-pairs for this purpose. For example, there is a construction related to the squared 
eigenfunctions of the Lax operator L [TU]-|12j. 

An alternative approach is based on the explicit formulas for the A- operators (see EH 
Ho]). Some version of this approach has been suggested in ^H]. in Section 3 we generalize the 
main idea of this work to the case of the Landau-Lifshitz equation. As a result, a deep relation- 
ship between the algebra A of quasilocal recursion operators and the algebra of multiplicators 
of the Lax structure becomes evident. In particular, the existence of two recursion operators 
related by elliptic curve equation follows from a similar property for multiplicators. 



1 The Krichever-Novikov equation. 

The Krichever-Novikov equation jHJ EI can be written in the form 

u h =u xxx -^ + ^, PW=0. (1.9) 

Denote by G\ the right hand side of (jl.9j) . The fifth order symmetry of (jl.9)l is given by 

u A u 2 5«§ 25 U2,u\ 45 u\ 5 u 3 25 u\ 5 u 2 5 P 2 5 „ 

G2 = W5 — 5 1 k n F — ~ H F — ^ F 5- H — U\F . 

u\ 2ui 2 u\ 8 uf 3 uf 6 3 u\ 18 u\ 9 
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The simplest three conserved densities of (jl.9|) are 

1 u\ 1 P 1 u\ 3 u\ 5 u\ IP 2 5 

= < g^l _ 19ttgtzj 7 u| 35 4 45 uf _ 259^ 35 
P3 «? 2 u\ 3 «f 9 uf 8 wf 36 18 u\ 

_14 «| 1 _ 14 W _ _7_P« _ 14 (7V) 
9 u? f 27n 6 27 uf 27 u? 9 r "I" 

In the paper [2] the forth order quasilocal recursion operator of the form 

Tlx = D 4 X + ai D 3 x + a 2 D 2 x + a 3 D x + a 4 + G^ 1 ^ + u,/^ 1 ^, 

ou ou 

was found. Here the coefficients a« are given by 

,u 2 a u\ u 3 4P 
ai = -4— , a 2 = 6— -2 - — , 



Mi Uf Ui 3 Mi 



«5 ^3,0^2 ^ 4P 2 4 7i| 10 8 U 2 

u\ u\ u\ u{ u\ 9 u\ 3 u\ 9 3 -uf 

The following statement can be verified directly. 

Theorem 1. There exists one more quasilocal recursion operator for ()1.9|) of the form 

U 2 = Dl + b x Dl + b 2 D x + 6 3 P' + hDl + b 5 D x + b 6 - \u x D~ l5 -^ 
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where 



Uo U3 P Un 

bi = -6— , b 2 = -9— -2^ + 21^, 
U\ Ui u{ u\ 

Ua U-iUi „Uo Uo P' 

b 3 = —11— + 60^ + 14P^ - 57^ - 3—, 

U\ U\ U'l U\[ Ui 

h = -4^ + 38^ + 22^1 + 99^ - 155^ + ^£ - 44P M ' 



U\ u\ u\ u\ u\ 3 u\ u\ 

3 u\ u\ 



b 5 = -2* + 29^? + 80P4 + ^P'^| - 104^1 - 70^1 + 241^% + 14^ 

Ml U\ U\ 3 ttj iq U\ U\ U\ 

20 «. _ 170 ^3 4 _ „| ^ _ 16 ^ _ «| 

3 u\ 3 3 m? u\ mi 3 

Mr UoUr 8 9 Mn u\u% u\ 142 Wo 28 , Mo UaU^Uo 

6 6 = - - 6^ + -P 2 ^ - 195^ + 6P^ + — P^ + -P'P^ + 101-4-^ 

Mi 9 u\ u\ 3 Ui 9 u\ u'l 

34 u A u 2 «g 28 38 W | 19 « 4 122 ^ nf nf 

+ y p yf ^"^ M2 + y p ^~y p ^~— p ^~ 10 ^ +22 ^ 

_ 178 p M3M| 14 (7V) 2 113 p/ M3M2 _ 2 p U5 _ 17 „U3 _ 4 p2 % _ 30 «4^ 
3 uf 9 1 3 wf 3 wf 3 «i 3 u\ u\ 

u 3 ui u 5 u 3 u h u\ IP' 2 8 P 3 4P"P 

+236^3 - 13^^ + 25^^ =- - — =-. 

uf u f 9 27 u? 9 u 2 

The operators 1Z\ and 1Z 2 are related by the following elliptic curve 

n 2 2 = n\- (pUx-e, (1.10) 



where 



i£ (^P»y _ 2 p'"p' + 2P^)p) ; 







16 
27 

243 



^ _ _(p") 3 _ _(p') 2 p( /y ) + pip" p>" _|_ 2p( /y ) P"P - P(P'") 2 j 



Remark 1. The relation (jl.lOj) is understood as an identity in the non-commutative field 
of pseudo-differential series [TS] of the form 



A = a m D™ + a m _iP™ 1 + ■ ■ • + a + a^ 1 D x 1 + a_ 2 D x 2 + ■ ■ ■ 
Here are local functions and the multiplication is defined by 

aD k x obD x n = a {bD^ +k + C l k D x {b)D k x +m - 1 + C 2 k D 2 x {b)D k+m - 2 + 



6 



where LmeZ and 



■ _ n(n - l)(n - 2) ■ ■ ■ (n - j + 1) 

Remark 2. It is easy to verify that <j) and 9 are constants for any polynomial P(u), where 
deg P < 4. Under Mobius transformations of the form 

au + (3 



u 



7« + 5 



in equation ()1.9|) the polynomial P(u) changes according to the same rule as in the differential 
uj = — ^^=. The expressions <p an d 9 are invariants with respect to the Mobius group action. 

v p ( u ) 

Remark 3. Of course, the ratio IZ3 = IZ^^ 1 satisfies equation ()0.4|) . It belongs to the 
noncommutative field of differential operator fractions JH]- Any element of this field can be 
written in the form P\P 2 l for some differential operators Pj. So, according to our definition, 
TZ 3 is a recursion operator of order 2. However, this operator is not quasilocal and it is unclear 
how to apply it even to the simplest commuting flow u x . 

The recursion operators presented above appear to be ratios 

1Z\ = T-LiHq 1 , 1Z 2 = 'H 2 H§ X 
of the following quasilocal Hamiltonian operators 

Hi = \{u 2 x Dl + /;>;) + {2u xxx u x - \u 2 xx - \P)D X + D x {2u xxx u x - \u\ x - \P) 
+G lJ D x T 1 G' 1 + u x D- l G 2 + G 2 D x l u x , 

H 2 = \tu\D\ + Dlul) + (3u xxx u x - fu 2 xx - P)Dl + P^a. - fu 2 xx - P) 
+hD x + D x h + G 1 D~ 1 G 2 + G 2 D~ 1 G 1 + u x D~ l G z + G 3 D~ 1 u x , 

where 

^ 2 ^ xxx / 2 \ xx ( 1 2 \ ^ 2 i lf 

h u xxxxx u x Qu xxxx u xx -\- — u xxx ~ (5.P 3Qu xx j -I- — (oP Qu xx j -\- — — -I- Uj-P , 

2 6 U x U x 6 U z x 

and G3 = 1Z\{Gi) = TZ 2 (u x ) is the seventh order symmetry of (jl.9|) : 



G 3 = m 7 - 7 7-9 ( 2P + 12w 3 wi - 27m 2 ) - h - — ~u 2 (2P - llu 2 ) 

U\ 6 u\ 2 u\ 2 u\ 



91 u 2 595 35«3 „, 2 1575,4 1813 k* 

"3 if 3 " 6 3 " is< (2F " P > ~ + "2T< 

.H2^ P < + «4( 6 pV - 5 P*) - I^(2P'V - 5PP') + If! 



7 P 7 7 P' 2 

- l -P"-+ l -P""u\-——. 
9 ui 9 18 Mi 
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2 The Landau-Lifshitz equation. 



In this section we consider the Landau-Lifshitz equation written in the form 

u t2 = -u xx + 2^ (ul - P(u)) + \P'{u) 
v t2 = Vxx + 2iP{vl-P{v))-\P'{v), 

where 

^ = (u-u) -1 , 

and P is an arbitrary fourth degree polynomial. The usual vectorial form of Landau-Lifshitz 
equation gives rise to a system of the form (|2.11|) after the stereographic projection (see Section 
3 for details). 

The third order symmetry of (|2.11)) is given by 

u t3 = u xxx - 6u x u xx ip + 6u 3 x ip 2 - \u x P"(v) - 3u x ipP'(v) - 6u x ip 2 P(v), 
v t3 = v xxx + 6v x v xx ip + Qvlip 2 - \v x P"{v) - 3v x ipP'(v) - 6v x ip 2 P(v). 

In case of the system of evolution equation (|2.11|) the symmetries and cosymmetries can be 
treated as two-dimensional vectors. We introduce the following notation for symmetries: 

G 1 = {u x ,v x )\ G 2 = (u t2 ,v t2 )\ G 3 = (u t3 ,v t3 y. (2.12) 

Thus Gij is j-th component of the symmetry Gi. The simplest cosymmetries are given by 

9i = ^ 2 (v x , -u x )\ g 2 = i?{v t2 , -u t2 )\ g 3 = ^ 2 (v t3 , -w t3 )*. (2.13) 

These cosymmetries are variational derivatives (J 2 , I 2 )* of the following conserved densities 

Pi = \ip{u x + v x ), p 2 = *P 2 (u x v x - P(v)) - \il)P'{v) - ±P"(v), 

P3 = \^ 2 {u x v xx - u xx v x ) + ij 3 u x v x (u x + v x ) - \u x ^(A^ 2 P(u) + P"{u) - 3ipP , (u)). 

The coefficients of operators F*,7l, and Ti. from f)0.4|) . f)0.6j) are matrices. For the quasilocal 
recursion operators the non-local terms can be written as AD~ l B t , where A and B are k x 2 
matrices whose columns are symmetries and cosymmetries, correspondingly. For quasilocal 
Hamiltonian operators the columns are symmetries for both A and B. 

Theorem 2. Equation (I2.11j) possesses the following quasilocal recursion operators: 



and 



*i = ° | - 2 ( G r 11 ° r 21 ) D? ( ^ 922 ) (2.14) 



Rll I V G 12 G 22 G 



r 32 



9n 9i2 



where G%j and gij are defined by (|2.12|) . ()2.13|) . and 

R l u = D 2 X - 4tfju x D x + 2^ 2 u x {v x + 3u x ) - 2^u xx - \P"{v) - AP{v)^ 2 - 2i/;P'(v), 
RI2 = D 2 X + Aij)v x D x + 2^ 2 v x {u x + 3v x ) + 2^v xx - \P"{v) - 4P(^ 2 - 2^P'{v), 
R 2 U =Dl- Qi/ju x D 2 x + (QiP 2 (3u 2 x - P(v)) - Qi/ju xx - \P"{v) - tyP'(v))D x + 4>{P"{v)u x - 2u xxx ) 
+4ip 3 u x (6P(v) - Qu 2 x -v 2 x - u x v x ) + ip 2 (2v x u xx - 2u x v xx + 9P'(v)u x + 18u x u xx ), 

R 2 2 = 2^ 2 (P(u) - u 2 x )D x - 4iP 3 (P(u)(u x - v x ) + u 2 x v x - u*) + ip 2 u x (P'(u) - 2u xx ), 
R 2 2l = -24, 2 {P{v) - v 2 x )D x + 4^(P(v)(u x - v x ) - v 2 x u x + v 3 x ) - ^ 2 v x (P'(v) - 2v xx ), 

Rh = -Dl - 6v x iPD 2 x + (Q^ 2 (P(v) - 3v 2 x ) - Qv xx i/j + 3^P'(v) + \P"(v))D x 

-4?p 3 v x (6v 2 x + v x u x + u 2 x ) + 2i) 2 (3v x P'(v) - 9v x v xx - u x v xx + v x u xx - \v x P'(u)) 
+\2^v x {P{v) + P{u)) + ip{v x P"{v) - 2v xxx ). 
Operators TZ± and 7I2 are related by the following elliptic curve equation 

n 2 2 -n 3 1 -ipTZ 1 -'dE = 0, (2.16) 
where E stands for the unity matrix, and 

if = I2ip 4 P(v)(P(u) - P(v)) + Ai) 3 {P'{v)P{u) - P{v)P'{u) - 3P(v)P'(v)) 
-i) 2 {2P"{v)P{v) + P'{v)P\u) + 3P'{v) 2 ) - i)P"{v)P'{v) - ±P"(v) 2 , 

$ = \Q^P(v)(P(v) - P{u)f + 8ip 5 (P(v) - P(u))(P(v)P'{u) + 3P'(v)P(v)) 

+^\P(v)P'{u) 2 - 3P'{v) 2 P{u) + 6P(v)(P'(u)P'(v) + 2P'{v) 2 ) + AP(v)P"(v)(P(v) - P{u))) 
+i) 2 P"{v){P'{v) 2 + \P"{v)P{v) + \P'{v)P'{u)) + \^P"{v) 2 P'{v) + j^P"{vf 
+^ 3 {P'{v) 2 P'{u) + lP"{v){P{v)P'{u) - P'{v)P{u)) + AP"{v)P{v)P'{v) + 2P'{v) 3 ). 

It is easy to verify that ip and t? are constants for any polynomial P, if degP < 4. They are 
invariants with respect to the Mobius group action just as in the case of the Krichever-Novikov 
equation (see Remark 2). 

Remark 4. Possibly the odd recursion operator found in [2H] is a rational function of 
recursion operators from Theorem 2. 

It turns out that the above recursion operators are ratios 
of the following quasilocal Hamiltonian operators 
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^12 | of ^11 ^21 A n _l / G*2l G22 



^ 7^21 J I G12 G22 J a V ^ n 



" 2 " ! ?4 7i| 2 J + 2 I G 12 G 22 G 32 ' D * I G21 G22 



where 



W12 = - 4^A* - 2— + 2v x + 61^ - - 2— - 4P(i;), 

1 1 9 «a!x 9 lP"(v) P'(v) 

H 2 U = 2(4 - P(u))D x + 2u,u„ - P'(u)u„ H% 2 = 2(t»J - P(v))D x + 2v,v xx - P»tJ, 



2 ip 2 ip J ip 

v 

_ 2 ^xx + _|_ 4. 4^1^ + + v x )(3P(v) - Au x v x ) + P'(v)(3u x + 6v x ), 

V 

+3w x P / (m) + - m / + + 8u x v xx + iv x u xx - QP'(v)u x 

ip ip 

+16u x v x (u x + v x )ip - \2ip(P(v)u x + P(u)v x ). 



3 Recursion operators and multiplicators. 

The original vector form of Landau-Lifshitz reads as follows 

Ut = Ux V xx + U x JU. (3.17) 

Here U = («i, w 2 , w 3 ), |U| = 1, symbol x stands for the vector product, and J = diag(p, q, r) 
is an arbitrary constant diagonal matrix. The usual way to represent equation (j3.17j) as a 
two-component system is to make use of the stereographic projection. The transformation 

U = (1 — uv, i + iuv, u + v) if), where i 2 = — 1, (3.18) 

coincides with the stereographic projection if we set v = —1/u. This transformation takes 
equation (|3.17|) into the system of the form (|2.11|) : 

u T = -u xx + 2^(u 2 x - P{u)) + \P'{u), 

(3.19) 



where t = —ir, 



v T = v xx + 2^{v 2 x -P{v))-\P'{v 
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Here and below the expression (A, B) denotes the standard Euclidean scalar product of the 
vectors A and B. Note that for equation (j3.19j) curve ()2.16|) has the form 



Recall the algebraic structure lying behind the elliptic Lax pair [21J for the Landau-Lifshitz 
equation. Let 

1\ /0 0\ / 1 0' 

e, (0 , e 2 = 1 , e 8 = -1 
-100/ \0 -1 0/ \ 

be the standard basis in the Lie algebra so(3). The Lax operator L for (|3.17|) is given by 

3 

L = D x - UiEu (3.20) 

3=1 

where 

1 / 7TT _ 1 



Ei = -e lV /l- P A 2 , E 2 = -e 2V /l-gA 2 , E 3 = - e 3 Vl - rA 2 . 
The operators Ai defining the Lax representations 

L u = [A, L) (3.21) 

for the commuting flows 

from the Landau-Lifshitz hierarchy belong to the Lie algebra Q generated by E 1; E 2 , E 3 . 
Lemma 1. The algebra Q is spanned by 

^Ex, ^E 2 , ^E 3 , ^E X> ^E 2) ^E 3 , i = 0,l,2,..., (3.22) 

where 

Ei = ^e lV / l-gA 2 v / T^ : A 2 , E 2 = ^e 2 v / l-j9A 2 v / T^7A 2 , E 3 = ^ e^l - pX^l - gA 2 

Our main observation is that the recursion operators are in one-to one correspondence with 
the multiplicators (see [2211221) of the algebra Q. 

Definition. A (scalar) function //(A) is called the multiplicator for the algebra Q if //(A) Q C 
Q. The order of pole of //(A) at A = is called the order of the multiplicator. 

It is easy to prove the following 



11 



Lemma 2. The set of all multiplicators for Q coincides with the polynomial ring generated 

by l, ^ 

^l(A) = "To"; /i 2 (A) - 



A 2' ^VV A s 

It is clear that 

^2 = (^i - - <?)Oi - 0- 

Thus the ring of multiplicators is isomorphic to the coordinate ring of an elliptic curve. 

The following construction establishes a correspondence between multiplicators and recur- 
sion operators. Let fi be a multiplicator of order k > 0. To find a relation between operators 
A n and A n+ k we use the following anzats 

A n +k = \xA n + R n} R n e G, ordR n < k. 

Substituting this into (j3.21|) and taking into account ()3.2U|) . we get 

±H^% = »±H<r>B 3 + [±u jEh R,] - (3.23) 

j = l j = l j = l 

Both sides of this relation belong to Q. Equating in ()3.23|) the coefficients of basis elements 



()3.22|1 . we find step by step unknown coefficients of R n and eventually an expression for #j n+fe) 



in terms of i.e. a recursion operator of order k. 
For the simplest multiplicator fii = A -2 , we have 

3 3 
3=1 3=1 

It is easy to verify that ()3.23|) is equivalent to the following identities: 

1. HW - F x U = 0, 

2. F I + UxM = 0, 

3. H( n+2 ) = Mj. + F x JU, 

where M = (M l5 M 2 , M 3 ), F = (Fi, F 2 , F 3 ). It follows from the first identity and the condition 
|U| = 1 that 

F = U x H (n) + / U. (3.24) 

To find / we note that the second identity implies (U, F x ) = 0. Substituting the expression 
(I3.24|) to this relation, we get 

f = D- 1 (v,n^x\j x ). 

The second identity can be rewritten as 

M = U x F x + mU = U(U, HjW) - + /UxU I + mU 
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for some function m. To find m, we substitute into the third identity the latter expression for 
M and take the scalar product by U of both sides of the relation thus obtained. Taking into 
account that (H (n+2) ,U) = 0, we get 

m = D^((JU,H^)-(\J x ,H x n) ))- 
Now the third identity produces the following recursion operator for equation (jH.17j) : 

H (n+2) = (U; H W )XJ _ H (n) + H W )U;c _ (U) Vx x H (n) )TJ x JJ x + J U)H W 

-(U x V xx + U x JV)D-\V, XJ X x HH) - V X D- 1 ((U x , H x n) ) - (JU, HW)) . 

(3.25) 

The recursion operator ()3.25|) . along with the bi-Hamiltonian structure of the Landau-Lifshitz 
equation, was originally discovered in [21|, whereas explicit formulas for higher symmetries were 
given earlier by Fuchssteiner [23] . 

For the second multiplicator /x 2 we set 

3 3 ^3 

3=1 j=l j=l 

Upon the substitution of (|3.26|) to (|3.23|) we get the following identities 

1. = F x U, 

2. F I = MxU, 

(3-27) 

3. M x - JHW = K x U, 

4. H( n+3 ) =K I + MxJU. 

The first two relations in (|3.27|) are analogous to (j2J), and therefore 

F = U x H (n) + fU, f = D-\V, H (n) x U x ), 

M = (U, H^ n) )U - H< n > + fVxV x + mU. 
It follows from (|3~2Tlh that 

K = Ux(M I - JH< n >) + k\J = (U, Hi n) )U x \J X - U x JH^ - U x 
-U a (U, HW x V x ) - ((U x , U X )V + V xx )f + mU x U x + kU. 

Functions m and /t can be found from the conditions 

(u, m x - jh^) = o, (u, h("+ 3 )) = o, 
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which yield 



k = D- 1 ((U,U X x JH<")) - (U,H^ x U x ) - i(U,HW x U X )((U, JU) + (U X ,U X )) 
+(U,Hi n) x JU)) - ±((U X ,U X ) - (U, JUJjD^HU.HW x U x ). 



H('"+ 3 ) = HS x U + xU, + (H™, U)U xU,- (U, x U !e )U 

+ (Hi n) , U)U x JU + (Hi n) , U)U x U^ - Hi n) x JU + (U, H x n) x JU)U 
-(U, Hi n) x U X )U X - 2(U, H< n > x U X )(U X , U X )U - \J X x JH^ - U x JH 
+(U, U x x JHW)U + (JH(»), U)U x U, - (U, HW x U^jU, 



One can verify that under transformation (J3.18)) the operators defined by formulas ()3.25|) 
and (|3.28|) become \{p + q + r)E — TZ\ and IZ2 correspondingly. 

Acknowledgments. The authors are grateful to Professors I.M. Krichever and W.K. Schief 
for useful discussions. The research was partially supported by: RFBR grant 05-01-00775, NSh 
grants 1716.2003.1 and 2044.2003.2. 



Thus the second recursion operator reads as 



-2(U, HW x U X )U XX - (U x U xx + U x JU)^- 1 ((Hi n) , U B ) - (JH^, U)) 
- ((U xx + §U(U X , U x )) x - §U X (U, JU)) D~ l (V, HW x U x ) 

+U X L>- 1 ((U,U X x JHW) _ (U,H<£ } x U x ) - |(U,HW x U X )(U, JU) 



(3.28) 
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